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' pH ' ■ Abstract 

We show the potential of the LHC to detect a CP-even Higgs boson of the NMSSM, 
h\ or /12, decaying into two rather light CP-odd Higgs bosons, a±, by exploiting the 
production mode based on Higgs-strahlung off 6-quarks, i.e., the channel pp — > bbh\^- 
We also consider the case of /12 —> h\h\ decays. It is found that these decays have 

■ dominant BRs over large regions of the NMSSM parameter space where tan/3 is large, 
, a condition which also favours the pp — >• 66/11,2 production rates. Further decays of 
' the light Higgs boson pairs (a\ and h\) into photon, muon, tau and 6 final states are 

■ also considered. The overall production and decay rates for these processes at inclusive 
level are sizable and should help extracting at least one Higgs boson signal over the 

\Q ' NMSSM parameter space. 

o ' 

1 Introduction 

^ ! The mechanism responsible for Electro- Weak Symmetry Breaking (EWSB) is still unknown. 

In the Standard Model (SM) and its extensions based on Supersymmetry (SUSY), such 
as the Minimal Supersymmetric Standard Model (MSSM) and the Next-to-Minimal Super- 
symmetric Standard Model (NMSSM), the Higgs mechanism is introduced for this primary 
purpose. Such a mechanism predicts the existence of at least one physical Higgs boson, 
which is a spin zero particle emerging from EWSB. While only one Higgs boson exists in 
the SM and five Higgs bosons in the MSSM, there are seven Higgs bosons in the NMSSM: 
three CP-even Higgses ^1,2,3 ( m h! < m h 2 < m h 3 ), two CP-odd Higgses a 12 (m ai < m a2 ) and 
two charged Higgses [T]. So, the latter has a phenomenologically richer Higgs sector than 
the two former scenarios. 

The NMSSM has two additional merits over the MSSM. On the one hand, it can solve 
the so-called ^-problem of the MSSM [2] in a natural way by introducing a new gauge singlet 
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field [3]. On the other hand, it can relieve the little hierarchy problem jU [5] since a SM- 
like scalar Higgs boson with mass less than the SM-like Higgs mass LEP limit is still quite 
naturally possible over some regions of the NMSSM parameter space. In fact, currently, the 
NMSSM can also explain a possible LEP excess and is definitely preferred by EW global 
fits. This happens when a SM-like Higgs boson of the NMSSM can unconventionally decay 
into two ai's with m ai < 2m b [6] (yet notice that this mass region is highly constrained by 
ALEPH [7\ and BaBar [8] data). In fact, there is also another possibility in the NMSSM, due 
to the fact that BR(ai — >■ 77) can be dominant and, as a result, BR(ai — > bb) is suppressed 
even though m ai > 2m b HD]- Finally, a CP-even Higgs state (hi or h 2 ) can naturally have 
a reduced couplings to the Z boson due to the mixing between the singlet and doublet Higgs 
fields, making it a natural possibility that the CP-even Higgs state of the NMSSM could 
have a mass less than the LEP limit on a SM-like Higgs mass. 

Probing the Higgs sector of the NMSSM is experimentally challenging, as it is not certain 
that we can always detect its physical states. There has been some work dedicated to explore 
the detectability of at least one Higgs boson of the NMSSM at the Large Hadron Collider 



and the Tevatron. In particular, some efforts have been made to extend the so-called 



'no-lose theorem' of the MSSM - stating that at least one Higgs boson of the MSSM will be 
discovered at the LHC via the usual SM-like production and decay channels throughout the 
entire MSSM parameter space [H] - to the case of the NMSSM P HOI H21 H31 HH HSl HE] . By 
assuming that Higgs-to-Higgs decays are kinematically not allowed, it was realised that at 
least one Higgs boson of the NMSSM will be discovered at the LHC. However, this theorem 
could be violated if Higgs-to-Higgs and/ or Higgs-to-SUSY particle (e.g., into neutralino pairs, 
yielding invisible Higgs signals) decays are kinematically accessible [TTt US]. 

So far, there is no conclusive evidence that the no-lose theorem can be confirmed in the 
context of the NMSSM. In order to establish the theorem for this SUSY scenario, Higgs-to- 
Higgs decays should definitely be taken into account though, in particular hi — > aia x . Such a 
decay can in fact be dominant in large regions of the NMSSM parameter space, for instance, 
for small Ak [9], and may not give Higgs signals with sufficient significance at the LHC. The 
importance of Higgs-to-Higgs decays in the context of the NMSSM has been emphasised over 
the years in much literature in all above respects, see, e.g., Refs. [51 ESI 1201 121] • Eventually, 
it was realised that Vector Boson Fusion (VBFjj could be a viable production channel to 
detect hi y 2 — > aiOi at the LHC, in which the Higgs pair decays into jjr + r~ [121 H21 122] . 
Some scope could also be afforded by a 4r signature in both VBF and Higgs-strahlung (off 
gauge bosons) [23]. The gluon-fusion channel too could be a means of accessing hi — » ai<2i 
decays, so long that the light CP-odd Higgs states both decay into four muons [21] or two 
muons and two taus [25] ■ Such results were all supported by simulations based on parton 
shower Monte Carlo (MC) programs and some level of detector response. 

In this paper, we want to investigate whether or not the no-lose theorem of the NMSSM 
at the LHC can possibly be reinforced by considering a Higgs production channel so far 
neglected, i.e., Higgs boson production in association with 6-quark pairs (aka Higgs-strahlung 
off 6-quark pairs). Notice that the twin process in which 6-quarks are replaced by t-quarks 

1 Hereafter, we consider 14 TeV as LHC energy. 
2 Which is dominated by W^+VF^-fusion over ZZ-one. 
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was discussed in [15], where it was found to be very subleading over the NMSSM parameter 
space. We will be looking at inclusive event rates in presence of various Higgs-to-Higgs 
decays, hi i2 — > a\ai and h 2 — > hihi, for hi and h 2 produced in association with fo-quark 
pairs. Notice that this production mode is in general the largest one in the NMSSM at large 
values of tan/3. We will also be studying the decay patterns of the lightest Higgs boson pairs, 
01O1 or h\h\, into different types of decay modes. 

This paper is organised as follows. In Sec. 2, we describe the parameter space scan 
performed. Inclusive event rates for the signals are explained in Sec. 3. Sec. 4 discusses 
possible signatures. Finally, we summarise and conclude in Sec. 5. 



2 Parameter Space Scan 

Due to the large number of parameters in the NMSSM, it is practically not feasible to do a 
comprehensive scan over all of them. Their number can however be reduced significantly by 
assuming certain conditions of unification. Since the mechanism of SUSY breaking is still un- 
known, to explore the NMSSM Higgs sector, we have performed a general scan in parameter 
space by fixing the soft SUSY breaking terms at high scale to reduce their contributions to 
the outputs of the parameter scans. Consequently, we are left with six independent inputs. 
Our parameter space is in particular defined through the Yukawa couplings A and n, the soft 
trilinear terms Ax and A K plus tan/3 (the ratio of the Vacuum Expectation Values (VEVs) 
of the two Higgs doublets) and /x cff — X(S) (where (S) is the vacuum expectation value of 
the Higgs singlet). 

We used here the fortran package NMSSMTools developed in Refs. [26j EZ| . This package 
computes the masses, couplings and decay widths of all the Higgs bosons of the NMSSM, 
including radiative corrections, in terms of its parameters at the EW scale. NMSSMTools 
also takes into account theoretical as well as experimental constraints from negative Higgs 
searches at LEP and the Tevatron as well as other contexts (S-physics, low energy 
experiments, etc.), including the unconventional channels relevant for the NMSSM. 

We have used the NMHDECAY code to scan over the six tree level parameters of the 
NMSSM Higgs sector in the following intervals: 

A : 0.0001 - 0.7, k: 0-0.65, tan/3 : 1.6 - 54, 
fi: 100 - 1000 GeV, A x : -1000 - +1000 GeV, A K :-10 - 0. 

Remaining soft terms, contributing at higher order level, which are fixed in the scan include: 

• m Q 3 = m U3 = m D:i = m Ls = m Es = 1 TeV, 

• A Ua = A Da = A E , = 1.2 TeV, 

• m,Q = mu = m D = m L = m E = 1 TeV, 

• Mi = M 2 = M 3 = 1.5 TeV. 

Notice that the sfermion mass parameters and the SU(2) gaugino mass parameter, M 2 , play 
crucial roles in constraining tan/3. Decreasing values of those parameters allow smaller values 
of tan/3 to pass experimental and theoretical constraints. In fact, when tan/3 is large the 
sfermion masses should be large to avoid the constraints coming from the muon anomalous 
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magnetic moment [29]. The dominant supersymmetric contribution at large tan/3 is due to 
chargino-sneutrino loop diagram [30]. Also, notice that the chargino masses depend strongly 
on M 2 . As mentioned above, we fixed the gaugino mass parameters and other SUSY breaking 
terms at high scale to reduce their contributions to the outputs of the parameter scans. 

In line with the assumptions made in [TJl H31 EH], the allowed decay modes for neutral 
NMSSM Higgs bosons ardj: 

h,a^-gg, h, a — >■ fi + , h,a— >t + t~, h, a — > bb, h,a—tti, 
h,a^-ss, h,a^>-cc, h — » W + W~ , h^-ZZ, 
h,a —> 77, h,a —> Zj, h, a — > Higgses, h,a — >■ sparticles. 

We have performed a random scan over 20 million points in the specified parameter space 
and required that m h2 < 300 GeV. The output of the scan, as stated earlier, contains masses, 
Branching Ratios (BRs) and couplings of the NMSSM Higgses, for all the successful points, 
which have passed the various experimental and theoretical constraints. 



3 Inclusive Event Rates 

For successful data points, we used CalcHEP [31] to determine the cross-sections for NMSSM 
Higgs production^. As the SUSY mass scales have been arbitrarily set well above the EW one 
(see above), the NMSSM Higgs production modes exploitable in simulations at the LHC are 
those involving couplings to heavy ordinary matter only. Amongst the production channels 
onset by the latter, we focus here on the processes 

gg, qq ->■ bb h x and gg, qq -> bb h 2 , (1) 

i.e., Higgs production in association with a 6-quark pair. This production mode is dominant 
at large tan j3. 

To a good approximation, at large tan /3, the tree level lightest neutral Higgs boson masses 
are given by the following expressions [33J: 
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3 Here, we use the label h{a) to signify any of the neutral CP-even(odd) Higgs bosons of the NMSSM. 

4 We adopt herein CTEQ6L [32] as parton distribution functions, with scale Q — y/§, the centre-of-mass 
energy at parton level, for all processes computed. Further, we have taken m^rrib) — 4.214 GeV for the 
(running) bottom-quark mass. 
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where 

(We reproduced here these tree level formulae mainly for guidance in interpreting the upcom- 
ing figures, the reader should recall though that NMSSMTools includes radiative corrections 
as well.) 

To probe the NMSSM parameter space, we have computed m hl and m^ 2 against each 
of the six parameters of the NMSSM (Figs. 1 and 2). As it is clear from the two figures, 
in our chosen parameter space regions, small values of A, k and fx e s are favoured whereas 
large values of tan/3 and positive values of A\ are the most compatible with theoretical 
and experimental data. The distribution over A K is uniform primarily because only small 
negative values of n are scanned over. 

Fig. 3 shows the correlations between all three Higgs masses, m ai , and rrih 2 . Since 
the successful points emerging from the scan have small values of A, k and also A K , only 
rather small values of m ai are allowed. It is remarkable that the smaller m ai the smaller 
and (two top-panes). In the bottom-pane of the same figure, for around 120 
GeV, m/jj can have values from just above up to slightly less than 120 GeV, showing the 
possibility that the two Higgs states can simultaneously have the same mass, m/ ll ~ m^ 2 . 
Notice also that the majority of points have between 115 GeV and 120 GeV, i.e., just 
above the LEP limit on a SM-like Higgs mass. 

The production times decay rates of hi and hi, in which h\ decays into two lighter a±s and 
/12 decays into either a pair of ai's or a pair of h^s, as functions of the Higgs masses and 
nih 2 (left-panes), tan/3 (middle-panes) and of the corresponding Higgs-to-Higgs decays BRs 
(right-panes), are shown in Fig. 4. In our choice of parameter space which has large tan/3 
we have noticed that the production rate of h\ in association with a bottom-antibottom, 
a(pp — > bbhi), is nearly constant, does not depend on the tree level parameters, while the 
production rate of h 2 , a(pp — > &6/12), is strongly dependent on tan/3 and other tree level 
parameters. In fact, notice that in the figure we multiply the production rates by the decay 
rates of Higgs-to-Higgs particles, which play crucial roles in the changes of the inclusive cross 
section. The two bottom middle-panes of the figure make clear that while large tan j3 values 
are a necessary condition for large production times decay rates of h 2 they are not a sufficient 
one, as most of the points accumulate at intermediate event rates. 

It is clear that Higgs-to-Higgs decays are dominant over a large area of NMSSM parameter 
space if Higgs-to-Higgs decays are kinematically allowed and so these decays should be taken 
seriously before claiming any validity (or otherwise) of the no-lose theorem for the NMSSM, 
see right-panes of Fig. 4. Fortunately, for considerable regions of parameter space, with 
different masses of h\ and h 2 , these production rates are sizable (up to 1000 fb or so), except 
for the case of h 2 — > h x h x where only few points have large production rates, due to smallness 
of BR(/i 2 — > h\hi) compared with BR(/ii j2 — > a\ai) in general. 

Fig. 5 displays the correlations between the three discussed production and decay pro- 
cesses. It is quite remarkable that the overall trend, despite an obvious spread also in the 
horizontal and vertical directions, is such that when one channel grows in event yield there 
is also another one which also does, hence opening up the possibility of the simultaneous 
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discovery of several Higgs states of the NMSSM (other than a\ also hi and h 2 ), an exciting 
prospect in order to distinguish the NMSSM Higgs sector from the MSSM one. 

In Fig. 6 we have calculated the signal rates of hi and hi through their cascade decays 
that finish with a x — > bb and/or a% — > t + t~ . It is shown that the signal rates are quite large, 
topping 1000 fb for hi and 100 fb for h 2 in case of 4b and 4r final states due to the fact that 
BR(ai — > bb) is dominant when m ai > 10 GeV and BR(ai — > t + t~) is dominant for m ai < 
10 GeV. The 2b plus 2r rates are one order of magnitude less than the former two due to 
the fact that only the parameter space points with m ai > 10 GeV have these final states in 
which BR(ai — > t + t~) is only about 10% of BR(ai — > bb). Overall, there are some regions 
of parameter space which have considerable signal rates that could be sufficient to discover 
the hi and h 2 through their aiai cascade decays at the LHC. 

The h 2 cascade decays ending with hi — > bb and/or hi — > t + t~ have less cross section 
(see Fig. 7). Only for nih 2 around 120 GeV the rates are quite sizable, topping 50 fb, 5 fb 
and 0.5 fb level for 4b, 2b plus 2r and 4r final states, respectively. 

As explained in Ref. [9], the BR(ai — > 77) can be dominant over a sizable region of 
NMSSM parameter space. This very peculiar phenomenon appears in this SUSY scenario 
(unlike the MSSM) because of the fact that a rather light CP-odd Higgs state can have a 
predominant singlet component and a very weak doublet one. As a consequence, all SM-like 
partial decay widths are heavily suppressed as they employ only the doublet component, 
except one: T(ai — > 77). This comes from the fact that the cii — > x + X~ coupling is not 
suppressed, as it is generated through the XHiH 2 S Lagrangian term and therefore implies 
no small mixing. Although the direct decay ai — > x + X~ is forbidden, the aforementioned 
coupling participates in the 0177 effective coupling. In hence, when BR(ai — > 77) is very 
large, no other SM-like BR can be. Hence, it makes sense to look at the scope of ai<ii — > 7777 
decays. The corresponding inclusive rates are found in Fig. 8 as functions of m^i and mi l2 
(two top-panes) for both hi — > aiai — > 47 and h 2 —> ai<ii — >• 47. Despite inclusive rates are 
never very large, it should be noticed a consistent population of points in the former case at 
mh x ^ 115 GeV yielding up to (9(1 fb) rates, with also a possibility of rates reaching up to 
100 fb for smaller m^, and in the latter case well spread out in m^ 2 from 115 to 300 GeV 
yielding some points between 0.1 and 1 fb. Moreover, we have shown in the same figure the 
inclusive results for hi — > aiai — > t + t~/j + u.~ and h 2 — > aiai — > t + t~ . The rates for hi 
can reach 1 fb -1 for various ranges of and roughly 0.5 fb -1 in case of h 2 for essentially 
any m h2 . 

Finally, notice that the cases hi — > ciidi — > [i + pT /i+ \i~ and h 2 — > aiai — > u. + /j 1 ~u, + fi~ 
contribute below the 0.01 fb level over the entire NMSSM parameter space, so we do not 
show the corresponding plots. 



4 Possible Signatures 

The production times decay rates presented in the previous section are inclusive results, 
whereby no cuts have been enforced on the final state particle^. Clearly, in order to detect 

5 Recall that we use a finite 6-quark mass, see Footnote 4. 
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the latter, a typical finite volume of an LHC detector has to be emulated. Further, in 
order to assess the significance of the signal yield, a background simulation (within the 
same detector region) has to eventually be carried out. Here, in the spirit of Ref. [15] , we 
would like to discuss the possible scope of the possible aforementioned signatures, without 
however venturing in such a complicated simulations. The key issue to be accessed is clearly 
whether one or more of the 6-quarks produced in association with the Higgs state h\ or 
/12 (henceforth called 'prompt' 6-quarks) in process Q ought to be tagged as such. The 
relevance of this should be clear from inspecting Fig. 9. The 6-quarks in the final state 
often emerges from the splitting of a gluon inside the proton, hence they can be at very low 
transverse momentum (denoted here by pr b )- To enforce vertex tagging with good efficiency, 
say Sb = 60%, a minimum px b value is always required, the lowest reasonable threshold 
being 15 GeV or so [34"t 135] . For the case of a single 6-tag the overall efficiency at very low 
Higgs masses (irrespectively of considering either hi or h 2 being produced) is some 2-3%, 
eventually growing to 14-15% for very massive objects. For the case of a double 6-tag, we 
are instead speaking or corresponding rates at the 1% to 8% level, respectively. Clearly then, 
the scope of the production and decay channels investigated in the previous section much 
depends on the Higgs mass produced and the decay patterns pursued. 

We reckon that for the 47 signature it should not be necessary to tag any of the 'prompt' 
6-quarks at all, as any of the (typically high transverse momentum and isolated) decay 
products of the a^s could act as trigger and the SM backgrounds (which would generally 
be induced by non-QCD processes) should not be prohibitively largqj. Regarding signatures 
with t's, for 4r and 2r2b one could certainly exploit a r trigger (both leptonic and hadronic) 
[341 135] . however (especially in case of hadronic r decays), it may be necessary to tag at 
least one 'prompt' 6-quark to suppress QCD backgrounds mimicking r — > hadrons. The 
case 2r2/i would clearly exploit a muon trigger instead. Finally, the case of a 46 signature 
of hi t 2 —> flidi and h 2 — >■ hxhi decays is totally unexplored, especially considering the fact 
that the entire final state would be made up of six quarks, i.e., with an unavoidable huge 
combinatorics and burdened by an extremely large pure QCD background. 

In essence, only a dedicated kinematical analysis of the decay products could in the end 
ascertain the true selection efficiency of a signature and its scope. What we can responsibly 
do here is to highlight three possible scenarios. Firstly, one whereby the signal rates in the 
proceeding section will not be reduced substantially after enforcing acceptance cuts: this is 
certainly applicable to 47 events emerging from a\ states with masses between 50 and 100 
GeV (where the BR(ai — > 77) is maximal, see bottom- left pane of Fig. 2 in [9]) and 4r 
events (with the heavy leptons decaying leptonically to electron and muons, which however 
induce a 1% suppression because of the consequent BRs). Secondly, one whereby all decay 
signatures involving (one or more) hadronic t's and 6's are reduced by a factor between 
7 and 50, depending on the produced Higgs mass, assuming a single tag only of 'prompt' 
6's. Thirdly, one whereby most possibly the 6b final state requires a double tag of 'prompt' 
6-quarks, reducing the signal yield by a factor between 20 and 100, depending on the h 12 

6 This is in fact very important in view of the fact that the 47 decay rate is the smallest one amongst 
those studied here. 
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Conclusions 

Searching for NMSSM Higgs states at the LHC is very complicated compared to the MSSM 
ones due to the dominance of Higgs-to-Higgs decays in large parameter space regions of 
the next-to-minimal SUSY model. This is the main reason why a no-lose theorem has not 
been confirmed yet in the context of the NMSSM. In view of this and following on previous 
work, where the case of VBF and Higgs-strahlung of W, Z bosons and t quarks was studied 
[T5] . we have found here that, at large values of tan/3, hi and h 2 production in association 
with bottom-antibottom pairs and decaying into lighter Higgses can have sizable signal 
rates in some regions of NMSSM parameter space, in a variety of decay patterns including 
photons, muons, tauons and 6-quarks themselves. We have verified this at the inclusive level 
and discussed what could happen in presence of acceptance cuts and consequent detector 
efficiencies. 

Clearly, in the end, only a dedicated decay analysis, in presence of not only acceptance 
but also selection cuts (the latter driven by the also necessary background assessment), will 
decree whether signal extraction is possible and through which signatures. However, our 
present study, alongside the findings of [T5], should eventually direct the NMSSM parameter 
space exploration where discovery significances can be found. In all circumstances, just like 
with other previous attempts at extracting NMSSM Higgs-to-Higgs signatures, evidence of 
those investigated here will require a rather large LHC luminosity sample, of (9(100 fb^ 1 ) or 
more. 
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7 Notice that for a\ masses comparable to typical transverse momentum thresholds of the decay products 
further severe reductions could occur, however, there is plenty of NMSSM parameter space giving sizable 
signals for heavier a\ states for all signatures considered here. 
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gure 1: The lightest two scalar Higgs masses and rrih 2 as functions of A, k and tan/3 
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gure 2: The lightest two scalar Higgs masses and nih 2 as functions of // e ff , A\ and A K . 
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Figure 3: The correlations between the lightest CP-odd Higgs mass, m ai and the lightest 
two CP-even Higgs masses, and nih 2 and between the latter two. 
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Figure 4: The rates for a(pp — > bbhi) BR(/ii — » aiai), cr(pp — >■ 66/12) BR(/i2 — >■ aiOi) and 
cr(pp — >■ 66/12) BR(/i 2 —> h\hi) as functions of corresponding Higgs masses, of tan/3 and of 
corresponding BRs. 
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Figure 5: The rates for a(pp — > bbhi) BR(/ii — > aia\) versus a(pp — > bbh 2 ) BR(h 2 — > aiai), 
for a(pp — > bbhi) BK(hi a\a\) versus a(pp bbh^) BR(/^2 — > h\h\) and for a(pp —> 
bbhi) BR(/i2 — > didi) versus a(pp —> 66^2) BR(/i2 — > h\h\). 
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Figure 6: The signal rates for a(pp — > bbhi) BR(/ii — > a\ai) and a(pp — > BR(/i2 — > 

a,ia\) times BR(aiai — > bbbb), BR(aiai — > bbr + T~) and BR(aiai — > t + t~t + t~) as functions 
of m/ ll and m h2 . 
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Figure 7: The signal rates for a(pp — > bbh^) BR(/i2 — > h\hi) times BR(/ii/ii — > bbbb), 
BR(hihi —> bbr + T~) and BR(hxhi — > t + t~t + t~) as functions of nih 2 . 
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Figure 8: The signal rates for a{pp — > bbh\) BR(/ii — > a\ai) and a(pp — > bbh 2 ) BR(/i 2 
aiai) times (BR(ai — > 77)) 2 and times BR(ai — > t + t~) BR(ai — > as functions of nih 

and m h2 . 
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Figure 9: The efficiency to tag one or two 'prompt' 5-quarks in the final state, given as 
percent ratio of the production cross section for pp — > bb Higgs (where Higgs can equally 
refer to an hi or h 2 state) after a the cut p Tb > 15 GeV over the total one (also including the 
6-tagging performances, e b and e\, respectively), as function of the Higgs boson mass. The 
distributions have been produced at parton level by using CalcHEP. Herein we use e& = 60%. 
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